Sensors 2014, 14, 3267-3283; doi:10.3390/sl40203267 



OPEN ACCESS 



sensors 

ISSN 1424-8220 

www.mdpi.com/joumal/sensors 

Article 

Modelling of Impulsional pH Variations Using ChemFET-Based 
Micr ode vices: Application to Hydrogen Peroxide Detection 

Abdou Karim Diallo Lyes Djeghlaf Jerome Launay ^'^ and Pierre Temple-Boyer 

^ CNRS, LAAS, 7 avenue du colonel Roche, F-31400 Toulouse, France; 

E-Mails: dialloabdoukarim@yahoo.fr (A.K.D.); ldjeghla@laas.fr (L.D.); jlaunay@laas.fr (J.L.) 
University of Toulouse; UPS; LAAS; F-31400 Toulouse, France 

* Author to whom correspondence should be addressed; E-Mail: temple@laas.fr; 
Tel.: +33-561-336-954; Fax: +33-561-336-208. 

Received: 30 December 2013; in revised form: 7 February 2014 /Accepted: 12 February 2014 / 
Published: 19 February 2014 

Abstract: This work presents the modelling of impulsional pH variations in microvolume 
related to water-based electrolysis and hydrogen peroxide electrochemical oxidation using 
an Electrochemical Field Effect Transistor (ElecFET) microdevice. This ElecFET device 
consists of a pH-Chemical FET (pH-ChemFET) with an integrated microelectrode around 
the dielectric gate area in order to trigger electrochemical reactions. Combining 
oxidation/reduction reactions on the microelectrode, water self-ionization and diffusion 
properties of associated chemical species, the model shows that the sensor response 
depends on the main influential parameters such as: (i) polarization parameters on the 
microelectrode, i.e., voltage (Vp) and time (tp); (ii) distance between the gate sensitive area 
and the microelectrode (d); and (iii) hydrogen peroxide concentration ([H2O2]). The model 
developed can predict the ElecFET response behaviour and creates new opportunities for 
H202-based enzymatic detection of biomolecules. 
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1. Introduction 

Hydrogen peroxide (H2O2) is an important chemical species and its determination is of great 
interest for food processing, industrial, clinical and biochemical applications. In the past years, many 
methods for H2O2 detection were published in the literature. These different techniques involve 
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spectrophotometry [1], chemiluminescence [2], titration [3] and fluorescence [4]. Nevertheless, they 
present considerable drawbacks, can be expensive and lack appropriate simplicity when used in 
specific applications. Electrochemistry is also a very interesting method for the detection of H2O2, 
leading to the development of associated sensors in liquid [5-8] and gaseous phases [9-11]. Indeed, 
amperometric or potentiometric electrochemical detection principles are advantageous over other 
methods due to their high sensitivity, fast response, low cost, simple instrumentation, compatibility 
with microtechnologies and related potentialities for miniaturization. The amperometric technique 
monitors redox phenomena in a conductive liquid/solid interface and has led to the development of 
electrochemical electrodes and cells [12,13] while the potentiometric one involves the detection of 
charges trapped at an insulating liquid/solid interface and has been responsible for the development of 
ion sensitive electrodes (ISE) and chemical field effect transistors (ChemFET) [13,14]. Both 
detection/transduction principles induced successful alternatives for the liquid phase analysis due to/in 
spite of different advantages/drawbacks in terms of (bio)chemical species detection, technological 
integration and data treatment. As a result, each of them was largely confined to specific applications 
according to their characteristics. Finally, the combination of amperometric and potentiometric 
techniques has also been very promising in terms of detection in liquid phase [15-19]. Such 
combination is associated with the functional integration of an electrochemical microelectrode and a 
pH-ChemFET on a single chip, leading to the realization of a pH-ChemFET based coulometric 
sensor-actuator system [15] also known as electrochemical field effect transistor (ElecFET) [19]. 

In parallel with the technological development of pH-ChemFET based microdevices, modelling 
investigations were also conducted. Based on the site-binding theory [20,21], such studies led to the 
development of complete behavioural models using SPICE [22], MATHEMATICA [23] or 
VHDL-AMS [24] softwares. From a general point of view, these different models focus on the study 
of electrolyte/insulator/semiconductor capacitive structures while taking into account the FET 
electrical behaviour due to well-known current- voltages equations [24,25]. Beyond the pH 
measurement, it is necessary to consider other phenomena occurring in liquid phase such as diffusion, 
migration, acid/basic chemical reactions, enzymatic biochemical reactions, .... This was performed for 
EnFET microsensors [26,27] as well as for pH-ChemFET-based coulometric sensor-actuator 
system [28,29]. Nevertheless, such modelling efforts have to be carried out in order to completely 
understand the ElecFET detection/transduction principles. 

In this paper, we report modelling of ElecFET microdevices in the frame of water electrolysis and 
hydrogen peroxide oxidation using the MATLAB™ software, studying main influential parameters 
(polarization voltage and time as well as integration level) while dealing with impulsional pH 
variations in microvolume and H2O2 detection. 

2. Experimental Section 

The ElecFET microdevice was developed through functional integration of a metallic 
microelectrode with a pH-sensitive chemical field effect transistor (pH-ChemFET) on a silicon 
chip [19]. Thus, by combining at the microscale amperometric production of water-based ions in 
solution and pH potentiometric detection, unusual electrochemical detection properties were 
evidenced [15-19] and should be modelled accordingly. 
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Since the pH-ChemFET theoretical behaviour is well known [14,24,25], its modelling was 
considered globally by taking into account the threshold voltage variations with pH (see Section 2.3). 
Nevertheless, in order to completely understand the ElecFET detection/transduction principles, its 
technological processes and the associated planar structure should be described (Figure 1) [19,30]. At 
first, N-channel, field effect transistors (FET) were fabricated on silicon wafers using a standard P-well 
technology. Thus, heavily-doped, N^-type source (S) and drain (D) regions (thickness: ~2 [im) were 
fabricated on a P-type silicon substrate to form the FET channel (length Lq: ~5 |im). Then, the 
pH-ChemFET gate was created. Typically, a 50 nm-thin silicon oxide Si02 layer was thermally grown 
on the silicon substrate in order to obtain optimal gate dielectric properties and a 50 nm-thin Si3N4 
layer was then deposited to provide pH detection properties. Finally, the metallic microelectrode was 
integrated on the device, next to the Si02/Si3N4 pH-ChemFET sensitive gate (standard distances: 
30-300 |im). Platinum was chosen since it is as an excellent catalyst for hydrogen peroxide 
oxidation [5,6,31]. Practically, platinum metallization (thickness: -200 nm) was used with a tantalum 
underlayer (thickness: -50 nm) to assure adhesion on Si3N4 surface. Overall, platinum microelectrode 
(area: -0.3 mm^) and pH-ChemFET (sensitive gate area: -4 x 10~^ mm^) were placed sufficiently close 
to one another in order to easily perform liquid phase analysis while keeping an adequate distance to 
avoid electrical interferences. 

Concerning the operating principle, a typical electrical configuration was used [19]. On one hand, 
the pH-ChemFET electrical bias was achieved through a constant drain-source voltage Vds and 
drain-source current Ids- Gate-source voltage Vqs was measured due to a reference electrode biasing 
the electrolyte to the mass (Vq = Vre = Veiectroiyte = 0). On the other hand, a suitable electrical 
polarization (voltage Vp, time tp) was applied on the integrated platinum microelectrode and a standard 
platinum counter-electrode in the solution was used to complete the set-up (Vce ~ Veiectroiyte = 0). 

Figure 1. Cross-section of the ElecFET device describing the one-dimensional model in 
radial coordinate (scales are not respected for easiness of presentation). 
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The ElecFET modelling takes into account the different chemical and physical phenomena 
occurring in the environment around the ElecFET microdevice: electrochemical production/ 
consumption of water-based acid/basic species (hydronium HsO^ and hydroxide 0H~ ions) due to 
redox phenomena on the integrated microelectrode, water self-ionization and diffusion phenomena in 
water medium. We assumed that (i) the transport by electromigration of HsO^ and 0H~ ions is 
negligible and (ii) any similar phenomenon is induced by other ions in the electrolyte. Similarly, 
transport by convection was neglected. As a result, modelling procedure is described as follows: 

• Study of the electrochemical production of water-based ions (HsO^ and 0H~) on the 
integrated microelectrode, 

• Study of the influence of water self-ionization, 

• Study of diffusion phenomena in liquid phase, determination of water-based ion 
concentration profiles and pH distribution around the pH-ChemFET sensitive gate, 

• Final determination of the pH-ChemFET potentiometric response (by considering only the 
typical "Nemst law" equation). 



2.1. Modelling of the Electrochemical Production of Water-Based Ions 

The ElecFET detection principles were first applied to water electrolysis phenomena triggered by 
an electrical polarization (voltage Vp, time tp) on the integrated microelectrode: 

• H2O oxidation (Vp > Eo+): 6 H2O — > O2 + 4 HsO^ + 4 e" 

• H2O reduction (Vp < Eo-): 4 H2O + 4e" — > 2 H2 + 4 OH" 

where Eo+ and Eo- are equilibrium potentials of oxidation and reduction reactions for water electrolysis 
respectively (Eo+ ^ 1.2 V and Eo- ^ -0.8 V [32]). 

The choice of a constant voltage to bias the integrated microelectrode was not compulsory in the 
frame of the ElecFET modelling: a constant current bias could have also been used. It was done in 
order to fit with ElecFET experimental characterisations previously performed [19]. Thus, according to 
laws of electrochemistry applied to a standard redox reaction (ox + ne~ ^ red), the current I on the 
electrode is related to the electrode potential V through the following equation [32]: 



I = nFS 



kox[red]exp 



anF , 



-V 



-kj.e^[ox]exp 



pnF 
RT 



V 



(1) 



where F is the Faraday constant (F = 96,485 C/mol), S is the electrode surface, kox and kred are the 
standard rate constants of oxidation and reduction respectively, a and P are the anodic and cathodic 
transfert coefficients respectively (a + P = 1), R is the ideal gas constant (R = 8.32 J/(K-mol)) and T is 
the absolute temperature (K). 

If the electrode potential V is higher than the oxidation potential E+, the reduction current can be 
neglected and Equation (1) becomes: 



I = nFSkQx[red]exp 



anF 



RT 



= nFSk+[red]exp -^(V-E+) 



VRT 



anF 



(2) 



Sensors 2014, 14 



3271 



Thus, in the case of monovalent ions, considering an elementary volume (surface S, thickness dr) 
the current I is given by: 

d[ox] 



I = eNSdr- 



dt 



(3) 



where N is the Avogadro number (N = 6.02 x 10^^). 

From Equations (2) and (3), temporal variations of the oxidant species concentration is finally 
given by: 



\^ = n-^[red]exp -—(V -E+) 
dt dr V RT 



(4) 



In the same way, if the electrode potential V is lower than the reduction potential E-, the temporal 
variations of the reducer species concentration is given by: 

anF 



d[red] k_ 
— — = n— [redjexp 
dt dr 



k_ =kredexp 



vRT 

anF 
RT 



(E_-V) 



(5) 



By applying Equations (4) and (5) to water oxidation and reduction respectively while considering 
that a = p = 0.5 [32], the temporal variations of HsO^ and 0H~ ion concentrations at the microelectrode 
surface (r = 0, cf. Figure 1) are finally given by: 



Vp > Eo+: 



Vp < Eo-: 



a[H30+] 
dt 



V 



J 



r=0 



k I 2F 
= G, =4^exp — (Vp -Eo+) 



dt 



't=0 



k I 2F 
= G_=4-f^exp— (Eo_-Vp) 
dr vRT 



(6) 



(7) 



where ko+ and ko- are standard rate constants of oxidation and reduction respectively (estimated 
respectively at 3 x 10"^ m/s and 10"^ m/s using the Tafel experimental method [32], results not shown), 
F is the Faraday constant (F = 96,485 C/mol), R is the ideal gas constant (R = 8.32 J/(K-mol)) and T is 
the absolute temperature (K). 

In presence of hydrogen peroxide in solution (concentration [H2O2]), the associated oxidation 
reaction should also be taken into account on the integrated microelectrode: 

H2O2 oxidation (Vp > Er): H2O2 + 2 H2O ^02 + 2 HsO^ + 2 e 

where Ei+ is the equilibrium potential for the H2O2 oxidation (Ei+ ^ 0.3 V for pH = 7 [32]). 

As previously stated, according to Equation (4), temporal variations of [HsO^] concentration at the 
microelectrode surface (r = 0, cf. Figure 1) are given by: 



Vp>E 



a[H30+] 



dt 



^r=0 



= Gi+ =2[H202]-|^exp — (Vp -E,^) 
dr VRT 



(8) 



where ki+ is the standard rate constant for the H2O2 oxidation (estimated experimentally at 

g 

5x10 m/s from the current lo obtained for Vp = Ei+ in agreement with the Tafel method [32], Figure 2), 
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F is the Faraday constant (F = 96,485 C/mol), R is the ideal gas constant (R = 8.32 J/(K mol)) and T is 
the absolute temperature (K). 

In our model, these equations were finally used to determine the electrochemical production of 
hydronium or hydroxide ions on the integrated microelectrode and its impact on the associated 
concentrations [HsO^] and [0H~]. 

Figure 2. Definition of the standard rate constant ki+ for the H2O2 oxidation using the 
Tafel method (S = 0.3 mm^ [H2O2] = 200 mM). 
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2.2. Influence of Water Self-Ionization 



In parallel, since the electrical polarization on the integrated microelectrode is responsible for the 
mass production of HsO^ or 0H~ ion, water is locally set out of equilibrium from a chemical point of 
view. According to water self-ionization (2 H2O HsO^ + 0H~), the kinetics of the [HsO^] and [0H~] 
concentrations are identical: 



d[H30+] d[OH"] 



(9) 



dt dt 

Introducing a constant parameter c that depends of the initial difference between the [HsO^] and 
[OH ] concentrations, the mathematical integration of Equation (9) gives: 

[H30^] = [0H"] + c <^ c = [H30^]-[0H"] (10) 

Thus, considering that the return to equilibrium is associated with a steady-state regime [29], the 
final concentrations of hydronium and hydroxide ions [HsO^Jf and [OH ]f are characterized by the 
following system (where Kw is the ionic product of water H2O): 

[H30+]f[0H-]f =K^ 

[HsO+Jf =[0H~]f +c ^ c = [H30+]f -[0H~]f 
Solving this equation system yields: 



(11) 
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[H3O+]?-c[H3O+]f-K^=0 



(12) 



Considering the positive root of Equation (12), the [H30"^]f and [OH ]f values are finally given by: 

(13) 



[H30+]f = — L w 



[OH-]f = 



(14) 



Finally, the global variations of the [HsO^] and [OH ] concentrations related to the water 
self-ionization reaction are expressed by: 

V([H30+]-[OH-])^ + 4K^ - [H3O+] - [OR-] 



A[H30+] = A[OH~] = 



(15) 



Thus, in the developed model, any variations of [HsO^] and [0H~] concentrations were 
counterbalanced according to Equation (10) in order to take into account water self-ionization 
phenomena in solution. 

It should be mentioned that the water self-ionization modelling was considered only in the case of 
pure water and therefore without considering any buffer properties. Such assumption is required to 
understand ElecFET detection/transduction principles. In practice, each solution of interest will 
have to be separately studied. Nevertheless, as demonstrated for urea-EnFET in the frame of 
haemodialysis [26,27,33], buffer effects will influence the ElecFET detection properties but will not 
drastically limit its operating principle. 

2.3. Modelling of Diffusion Phenomena in Watery Phase 

The diffusion model was adapted from a previous one developed for the EnFET-based microsensors 
and based on a finite element model implemented in MATLAB™ software [26,27]. It is associated 
with Tick's diffusion equation, assuming a one-dimensional model in spherical coordinates (Figure 1). 
The variable parameter is therefore the radius of the sphere r: 



aC(r,t) _^ 1 d 
dt dr 



dr J 



+ G 



(16) 



where C(r,t) describes the concentration of HsO^ or 0H~ ion for the studied case, D is the associated 
diffusion coefficient and G represents the chemical species generation according to the associated 
electrochemical reaction. 

As a result, for the different studied cases. Equation (11) gives: 



H,0 oxidation (Vp > Eo.): Ml^JM =D^.\^ 



.2 g[H30+](r,t) 
dr 



(17) 



H.0 reduction (Vp<Eo.): ™^^=D^^_-L| 



2 a[0H-](r,t) 

^ dr 



+ G. 



(18) 
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H,0, oxidation (Vp > Ei.): MsO^lM = Yy^^,\^ 



'^2e[H30^](r,t)' 



dv 



+ Gi+ (19) 



In order to solve diffusion equations, the following initial and boundaries conditions were chosen: 

[H3O+](r,0) = 10"P"o ^ [OH"](r,0) = K^lO+P^o 
'aC(r,t)\ 



at 



'dC{xX)\ (20) 



Thus, initial conditions assume spatially uniform concentrations of water-based ions and boundary 
conditions assume that no flux goes through the physical barrier of the sensor surface (r = 0) or far 
from it (r = rmax = 3 cm). By solving the mass transport equations system (Equations (17), (18) and/or 
(19)) using MATLAB™ software, the different ion concentration distributions [H30^](r,t) and 
[OH~](r,t) were defined and the pH(r,t) function was deduced consequently. 

2.4. Modelling of the pH-ChemFET Electrical Behaviour 

Finally, concerning the pH-ChemFET electrical behaviour, modelling of the drain-source current 
Ids as a function of gate-source and drain-source voltages Vqs and Vds, pH and any other interfering 
parameters was not performed. Indeed, since the pH-ChemFET device is well known at a theoretical 
level and was thoroughly studied in previous works [14,24,25], it was decided to take only into 
account its global detection properties associated with the Nemst's law. Thus, since the channel length 
Lg (typical values around 5 |Lim [30]) is shorter than the distance d between the integrated 
microelectrode and the pH-ChemFET gate sensitive area (Figure 1), the pH-ChemFET threshold 
voltage Vt variation was estimated in a very simple way according to the following equation: 

5VT(t) = So(pH(d,t)-pHo) (21) 

where so is the pH-ChemFET sensitivity and d is the distance between the integrated microelectrode 
and the pH-ChemFET gate sensitive area (Figure 1). 

3. Results and Discussion 

Hereinafter, the ElecFET model is studied using the following parameters: 

• T = 300 K and So = 60 mV/pH 

• Kw= 10"^^(mol/L)2and pHo = 7 

• Dh3o+ = 9.3 X 10"^ cm^/s and Dqh- = 5.3 x 10"^ cmVs 

• S = 0.3 mm^ and d = 30 |im (if not indicated) 

3.1. Modelling of Water Electrolysis Phenomena 

The first study concerned modelling of a positive bias on the integrated microelectrode (polarization 
voltage Vp = 1.23 V and polarization time tp = 5 s). In this case, water oxidation occurs on the 
integrated microelectrode, triggering the local production of HsO^ ions since Vp is higher than the 
water oxidation equilibrium potential Eo+. Figure 3 show the temporal and spatial variations of 
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hydronium (HsO^) and hydroxide (OH ) ion concentrations. As expected, the [HsO^] concentration 
(respectively, the [OH ] concentration) increases (respectively, decreases) very quickly by more than 
four decades (from 10"^ mol/L to 7 x lO "^ mol/L) and gradually takes its initial value (10~^ mol/L) as 
soon as the electrical bias on the integrated microelectrode is turned off (Figure 3a,c). Furthermore, it 
is obvious that these variations occur at millimetric distance (Figure 3b,d). This demonstrates that 
localized impulsional pH variations can be effectively obtained near the integrated microelectrode 
surface and can be detected at the adjacent pH-sensitive surface. 



Figure 3. (a,c) Temporal (r = d = 30 |Lim) and (b^d) spatial variations (t = 12 s) of the 
[HsO^] and [0H~] concentrations for a positive bias on the integrated microelectrode 
(Vp= 1.23 Vandtp = 5 s). 
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In the case of negative bias (polarization voltage Vp = -0.85 V and polarization time tp = 5 s), water 
reduction is taking place on the integrated microelectrode, triggering the local production of 0H~ ions 
since the Vp is lower than the water reduction equilibrium potential Eq-. Figure 4 illustrate respectively 
the temporal and spatial evolutions of the 0H~ ion concentration. An [OH ] increase (from 10 ^ mol/L 
to 4 X 10^^ mol/L) is noticed, followed by a return to equilibrium (10^^ mol/L) when the electrical bias 
is turned off. Of course, the HsO^ ion concentration follows opposite variations in agreement with the 
water acid/basic equilibrium (figures not shown). Impulsional pH variations are therefore evidenced at 
the millimetric scale, confirming finally that varied localized phenomena can be effectively obtained 
according to the polarization sign. 
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Subsequently, we have studied the influences of the main parameters by focusing on the 
pH-ChemFET threshold voltage variations 5Vt. Figure 5 shows the associated temporal variations for 
different polarization voltages applied on the integrated microelectrode (Vp = 1.21, 1.23, 1.25, 1.27 
and 1.29 V) and a given polarization time tp = 5 s. 



Figure 4. (a) Temporal (r = d = 30 |Lim) and (b) spatial variations (t = 12 s) of the [OH ] 
concentration for a negative bias on the integrated microelectrode (Vp = -0.85 V and 
tp = 5 s). 
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Figure 5. Temporal variations of the pH-ChemFET threshold voltage for different 
polarization voltages on the integrated microelectrode (Vp = 1.21, 1.23, 1.25, 1.27 and 
1 .29 V) and a given polarization time (tp = 5 s). 
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As expected, the polarization voltage increase is responsible for a local pH decrease and 
consequently a pH-ChemFET threshold voltage decrease. This behaviour corroborates the 
Butler- Volmer theory adapted to water oxidation (Equation (1)). Finally, since HsO^ ion production 
kinetics is an exponential function of Vp, an increase in polarization voltage induces a quasi-linear 
decrease of the minimum voltage. 
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Then, the influence of polarization time tp was studied. Figure 6 illustrates the temporal variations 
of pH and associated pH-ChemFET threshold voltage for different polarization times (tp = 0.2, 1,7, 15 
and 20 s) while keeping constant the polarization voltage on the integrated microelectrode 
(Vp = 1.23 V). In agreement with the electrochemical theory (Equation (1)), the polarization time 
increase is responsible for a local pH decrease and therefore a pH-ChemFET threshold voltage 
decrease. Nevertheless, pH variations are lower and tend to reach saturation. 

Figure 6. Temporal variations of the pH-ChemFET threshold voltage for a given 
polarization voltage on the integrated microelectrode (Vp = 1.23 V) and different 
polarization times (tp = 0.2, 1,7, 15 and 20 s). 
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Figure 7. Temporal variations of the pH-ChemFET threshold voltage for different 
distances between the integrated microelectrode and the pH-sensitive gate (d = 30, 90, 150 
and 210 |Lim) and a given polarization step (Vp = 1.23 V and tp = 5 s). 
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Since diffusion phenomena control the ElecFET detection principles, spatial configuration should 
influence detection properties. To demonstrate such assumption. Figure 7 shows temporal variations of 
pH-ChemFET threshold voltage for different distances between the integrated microelectrode and the 
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pH-ChemFET gate sensitive area (d = 30, 90, 150 and 210 |im) and for a given microelectrode 
polarization (Vp = 1.23 V, tp = 5 s). It appears that the distance diminution is responsible for a local pH 
and therefore a pH-ChemFET threshold voltage decrease. This demonstrates that the ElecFET 
detection principles are effectively dependent on the spatial integration of both microelectrode and 
pH-ChemFET microdevices. For the highest integration levels (d « 50 |Lim), very important localized 
impulsional pH variations are obtained for a given microelectrode polarization step (Vp, tp). On the 
contrary, for the lowest integration levels (d » 200 |Lim), phenomena are strongly attenuated, leading 
to low temporal pH wave and, consequently low ElecFET response. Since such pH variations were 
found to occur at millimetric distance (cf. Figure 3b,d), ElecFET detection potentialities is no longer 
possible when d parameter is higher than 1,000 |im, i.e., when the functional integration between the 
microelectrode and the pH-ChemFET device is ineffective. 

3.2. Modelling of Hydrogen Peroxide Oxidation Phenomena 

Since hydrogen peroxide detection is of great interest and H2O2 was successfully studied using an 
ElecFET-based microdevice [17,19], our model was studied by taking into account the H2O2 oxidation 
on the integrated microelectrode. In this case, according to the hydrogen peroxide oxidation reaction, a 
positive bias on the integrated microelectrode is responsible for the electrochemical production of 
HsO^ ions and therefore a local pH decrease in solution (cf. Section 2.1). Since H2O and H2O2 
molecules have similar electrochemical behaviours, impulsional pH variations localized at the 
microscale (typical dimension less than 1 mm) were also obtained in H202-rich solutions as soon as an 
appropriate polarization step (Vp, tp) was applied (result not shown). In this case, phenomena occur at 
lower voltage according to the associated equilibrium potential Ei+ value (Ei+ ~ 0.3 V for pH = 7 [32]). 
This is clearly evidenced on Figure 8 illustrating temporal variations of pH-ChemFET threshold 
voltage for different polarization voltages (Vp = 0.32, 0.35, 0.38, 0.42 and 0.45 V) and different 
polarization times (tp = 1, 3, 5, 10 and 20 s), considering an H2O2 concentration of 45 mM. In this case, 
local impulsional pH variations are not so obvious. Indeed, during the polarization duration (tp = 30 s), 
hydrogen peroxide oxidation is limited by the H2O2 diffusion phenomena towards the microelectrode 
surface. As a result, pH variations are no longer impulsional and higher polarization durations are 
required to reach a steady-state regime depending on the polarization voltage Vp (as well as on the 
[H2O2] concentration, see hereafter). Then, as previously stated, since the H2O2 oxidation reaction is 
controlled by the Butler- Volmer theory (Equation (3)), the polarization voltage Vp increase is 
responsible for a quasi-linear threshold voltage decrease while the polarization time tp increase is 
associated to saturation phenomena. 

Finally, in the case of hydrogen peroxide detection, the Butler- Volmer theory emphasizes on the 
influence of the [H2O2] concentration on HsO^ ion production kinetics (Equation (3)). In order to study 
this influence. Figure 9 shows temporal variations of the pH-ChemFET threshold voltage for different 
H2O2 concentrations ([H2O2] = 1, 10, 30, 60, 100 and 150 mM) and a given polarization step 
(Vp = 0.35 V and tp = 30 s). As expected, it is obvious that the hydrogen peroxide concentration increase 
is responsible for a local pH decrease and therefore a pH-ChemFET threshold voltage decrease. 
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Figure 8. Temporal variations of the pH-ChemFET threshold voltage for different 
polarization steps, (a) Vp = 0.32, 0.35, 0.38, 0.42 and 0.45 V, (b) tp = 1, 3, 5, 10 and 20 s, 
and for a given H2O2 concentration ([H2O2] = 45 mM). 




Figure 9. Temporal variations of the pH-ChemFET threshold voltage for different H2O2 
concentrations ([H2O2] = 1, 10, 30, 60, 100 and 150 mM) and a given polarization step 
(Vp = 0.35 Vandtp = 30 s). 
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In order to quantify the pH variation, its minimal value and more precisely the associated minimal 
threshold voltage value have been studied according to [H2O2] concentration for two integration levels, 
i.e., for two distances (d = 30 and 210 |im) between the integrated microelectrode and the 
pH-ChemFET gate sensitive area (Figure 10 respectively). Simulation results exhibit linear variations 
for several concentration decades (sensitivity ~ 60 mV/decade). Such Nemstian sensitivity is related to 
linear variations of the [HsO^] concentration with the [H2O2] concentration at the microelectrode 
surface (cf. Equation (3)). This result demonstrates that the pH-ElecFET microdevice can be 
effectively used for the potentiometric detection of hydrogen peroxide. Furthermore, by increasing 
the polarization voltage Vp and/or the integration level, H2O2 detection range and detection limit can 
be improved. 
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Figure 10. H2O2 analytical responses of the ElecFET microdevice (Vp = 0.32, 0.35, 0.38, 
0.42, 0.45 V and tp = 30 s) for two different distances between the integrated 
microelectrode and the pH-sensitive gate: (a) d = 30 |xm and (b) d = 210 |xm. 
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4. Conclusions 

Using MATLAB™ software, we have investigated the modelling of the ElecFET microdevice in 
the case of water electrolysis and hydrogen peroxide oxidation, taking into account chemical and 
electrochemical and acid/basic reactions in water phase as well as diffusion phenomena of water-based 
ions. In the first case, results show that the ElecFET phenomena depend strongly on the polarization 
parameters on the microelectrode (voltage Vp and time tp) as well as on dimensional characteristics 
related to the ElecFET functional integration (distance between the microelectrode and the 
pH-ChemFET gate sensitive gate area d). On one hand, according to the polarization sign, oxidation 
and reduction reactions of water are responsible for impulsional pH variations localized at the 
microscale (typical dimension: 1 mm and less). On the other hand, amplified phenomena are obtained 
when the polarization conditions and/or the integration level are increased. 

In the second case, the ElecFET microdevice was shown to be compatible with hydrogen peroxide 
detection while using an indirect measurement technique based on the pH-related H2O2 oxidation 
products. Amplification phenomena were shown to have no influence on the detection sensitivity and 
Nemstian responses were finally evidenced for several concentration decades. Nevertheless, detection 
range variations and improved detection limits were still possible by increasing polarization conditions 
and/or integration level. 

It should be mentioned that the ElecFET modelling was performed for solutions based on pure 
water and therefore without considering any buffer properties and/or any interfering pH variations. 
This choice was required to ElecFET detection/transduction principles, but it prevents us from 
comparing fully modelling and experimental results. As far as real samples are concerned, buffer 
properties will have to be considered to determine pH distributions in the ElecFET environment. Thus, 
the modelling will have to be adapted to each different buffer solution of interest. This can induce 
contradictory phenomena for the ElecFET behaviour. On one side, buffer properties can impede in 
some extents local pH changes, deteriorating detection properties in terms of sensitivity and/or 
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detection range. On the other side, by impeding pH variations in solution, they can improve detection 
performances in terms of selectivity (especially to pH). 

Nevertheless, the ElecFET microdevice is finally very promising in terms of applications, and 
modelling enables a complete comprehension of its detection principle based on a combination of 
amperometry and potentiometry. It allows a combination of the pH-ChemFET-metry technique and 
redox phenomena. Further studies will be focused on the modelling of the detection of H202-related 
biomolecules such as glucose as well as lactate, urate and/or glutamate ions by taking into account 
associated enzymatic reactions in buffered solution. 
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